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Tine,  Frequency,  and 
Space  Domains 

Volume  Only  (Direct)  Acoustic 
Paths 
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^Acoustic  propagation  over  two  purely  refracted  undersea  paths  is 
characterized  by  coherence  lengths  in  the  time,  frequency,  and  space  domains. 
Intensity  fluctuations  are  quantified  by  the  variance  of  log  intensity. 
Measurements  have  been  carried  out  on  1.5  and  18  nmi  paths  at  the  Azores 
Fixed  Acoustic  Range  (AFAR)  over  the  decade  of  acoustic  frequencies  from  500 
to  5000  Hz.  The  18  nmi  path  breaks  up  into  "micropaths,)’  resulting  in  self¬ 
interference  of  the  deterministic  ray.  _ 
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VOLUME  PROPAGATION  EXPERIMENTS  AT  THE  AZORES  FIXED 

Acosirric  range 


INTRODUCTION 


The  Azores  Fixed  Acoustic  Range  (AFAR)  was  an  international  research 
and  development  project  in  underwater  acoustics  in  which  eight  member 
nations  of  NATO  participated.1 

The  range  consisted  of  a  transmitting/receiving  array  mounted 
JO  m  above  the  ocean  bottom  at  a  depth  of  54C  m  (figure  1),  two  similarly 
constructed  receiving  arrays  at  depths  of  256  and  756  m,  and  an  omni¬ 
directional  hydrophone  at  a  depth  of  525  m.  The  three  array  positions 
formed  an  almost  equilateral  triangle  18  nmi  on  a  side.  The  hydrophone  was 
located  1.5  nmi  from  the  source  array.  Each  array  was  located  on  a 
mountain  top  plateau  and  the  valley  depth  between  these  terminals  reached 
2500  m.  The  arrays  were  steerable  in  azimuth  and  elevation,  and  spanned  the 
decade  of  acoustic  frequencies  from  approximately  500  to  5000  Hz. 

(hiring  the  first  half  of  the  1970' s,  a  variety  of  acoustic  experiments 
were  performed  at  AFAR  to  measure  the  effects  of  the  ocean  surface,  bottom, 
and  volume  on  acoustic  propagation  in  a  bistatic  situation.  The  published 
results  of  these  experiments  are  summarized  in  "The  AFAR  Project  - 
A  Valediction."2 

This  report  deals  with  a  series  of  measurements  carried  out  over 
volume-only  (direct)  acoustic  paths.  The  ocean  volume  is  common  to  all 
paths  between  source  and  receiver  and  its  effect  on  propagation  must  be 
measured  and  understood  if  we  hope  to  understand  undersea  acoustic 
propagation  in  general. 

There  were  two  direct  acoustic  paths  available  at  AFAR  that  could  be 
isolated  in  the  time  domain  from  boundary  reflected  paths.  Figure  2  shows 
a  ray  tracing  of  the  18  nmi  direct  path  linking  the  source  array  and  the 
southernmost  receiving  array.  The  direct  ray  was  confined  to  the  range 
of  depths  from  400  to  800  m  and  to  ray  angles  (referenced  to  the 
horizontal)  under  2  deg.  The  existence  of  this  path  is  due  to  a  positive 
"hump"  in  the  vertical  sound  speed  profile,  centered  at  a  depth  of  approxi¬ 
mately  1000  m,  possibly  caused  by  an  intrusive  outflow  from  the  Mediterranean 
Sea.  The  second  direct  path  at  AFAR  was  between  the  source  array  and 
the  omnidirectional  hydrophone  1.5  nmi  away.  The  experimental  results 
presented  below  show  that  the  two  paths  and  the  available  acoustic  frequency 
range  allow  investigation  of  both  weak  and  strong  scatter  propagation 
conditions . 
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Figure  1.  Overall  Front  View  of  an  AFAR 
Antenna  Tower 
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Figure  2a.  Average  of  57  Sound  Speed  Profiles  Measured  at 
AFAR  During  March  1975  Experiments 
(upper  840  m) 
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Figure  2b.  AFAR  18-nai  Direct  Ray  From  Speed 
Profile  of  Figure  2a 


Figure  2.  Direct  Path  Ray  Trace  and 
Sound  Speed  Profile 


SIGNAL  GENERATION  AND  DATA  ACQUISITION 


The  signaling  techniques  and  hardware  utilized  at  AFAR  are  described 
in  detail  in  the  "Exploitation  of  the  Azores  Fixed  Acoustic  Range  (AFAR) 
through  May  1973."* 

Briefly,  baseband  samples  of  the  complex  envelopes  of  signals  to  be 
transmitted  wire  stored  in  a  computer  memory.  The  samples  *  ere*  read  from 
memory  at  selectable  rates,  quadrature  modulated  to  the  desired  center 
frequencies,  digitally  filtered,  summed,  converted  to  analog,  and  distri¬ 
buted  to  the  appropriate  power  amplifier/acoustic  antenna  for  transmission. 

Received  signals  were  digitized,  quadrature  demodulated  to  baseband, 
digitally  filtered,  down  sampled,  and  recorded  on  magnetic  tape.  The 
quadrature  demodulation  to  baseband  minimized  the  required  recording 

sample  rate  and  also  facilitated  post -experiment  data  reduction. 

The  transmitter  and  receiver  were  slaved  to  a  common  frequency 
standard . el iminating  timing  problems,  part icularly  in  phase  dependent 
computations. 

Reference  to  data  reduction  in  this  report  will  assume  that  all 
operations  are  performed  on  samples  of  the  complex  envelopes  of  received 
signals. 


THE  JOINT  OCF.ANOGRAPHIC  EXPERIMENT 


The  task  of  performing  acoustic  measurements  under  all  bottom,  sur¬ 
face,  and  volume  conditions,  as  well  as  transmitter  and  receiver 
positions  and  speeds,  is  clearly  beyond  reach.  It  is,  therefore,  neces¬ 
sary  to  understand  the  physics  of  propagation  and,  hopefully,  to  develop 
a  model  that  will  allow  the  prediction  of  propagation  characteristic*  for 
any  combination  of  environment  and  geometry.  The  controlled  geometry  at 
AFAR. as  well  as  the  existence  of  a  strongly  scattered  volume  only  path, 
has  provided  a  unique  set  of  measurements  for  testing  volume  propagation 
models. 

During  March  197S,  an  international  acoustic-oceanographic 
experiment  was  conducted  at  AFAR  to  obtain  simultaneous,  detailed 

environmental  and  acoustic  data,  and  subsequently  rredict  the 
propagation  conditions  from  the  environmental  conditions.  Individual 

results  presented  in  this  report  are  primarily  from  the  March  1975  opera¬ 
tion  and  are  representative  of  other  AFAR  experimental  results. 
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PROPAGATION  MODELS 


At  the  time  the  AFAR  measurements  were  Mile,  no  propagation  models 
(of  which  the  author  is  aware)  existed  for  strong  scattering  by  random 
Inhomogeneities  in  the  ocean  volume  that  simultaneously  incorporated  the 
effects  of  anisotropy,  statistical  inhomogeneity,  and  the  deterministic 
sound  channel.* 1  Recent  work,  primarily  by  a  subgroup  of  the  JASON 
Committee,*1  has  had  remarkable  success  in  predicting  the  measured 
acoustics  using  oceanographic  data  obtained  during  these  experiments. 

The  model  includes  the  effects  of  anistropy,  inhomogeneity,  and  a  sound 
channel,  as  well  as  internal  wave  spectra. 


CHARACTERIZATION  OF  THE  PROPAGATION  CHANNEL 


We  have  chosen  to  quantify  the  acoustic  propagation  characteristics 
of  the  ocean  volume  in  terms  of  the  second  order  statistics  of  its  re¬ 
sponse  to  a  sine  wave  excitation,  i.e.,  the  time  varying  transfer  func¬ 
tion  H(f, t,r, ) .  When  H(f,t,r)  is  wide  sense  stationary  in  the  three 
variables  of  t  (time),  f  (frequency),  and  r  (space),  we  can  write  the 
second  order  statistic  as 


R(t1't2*  fl*f2'  rl'r2}  ’  H(tl*fl»rl)  H#ft2*f2,T2)  *  (1) 

When  we  suppress  the  mean  values  in  comjxjtation,  we  will  refer  to  R 
as  a  covariance  function. 

The  assumption  of  stationarity  can  not  be  Justified  a  priori  for  the 
ocean  volume.  We,  therefore,  make  a  less  restrictive  assumption 

R  -  Rftj-tj.fj-fj.Tj-rj.t^fj.r^  ,  (2) 

where  R  is  assiaaed  stationary  over  some  narrow  region  about  the  epoch, 
tQ,  center  frequency,  fQ,  and  spatial  point,  rQ. 

Finally,  using  t,  n,  and  S  as  shift  variables,  we  can  write 

R(T,n,S;  tQ,  fQ,r0)  -  H(t0M/2.f0*r/2, 

r  ♦  5/JJ'  H~«(t  -T/i.f  -n/27*vs7?)  •  (3) 

O  o  o  o 
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While  the  correlation  function  R  (  )  is  useful  in  quantifying  a  prop* 
gation  channel's  effect  on  signals,  the  physics  of  the  channel  are  most 
easily  understood  in  terms  of  Fourier  transforms  of  R  on  each  shift  vari¬ 
able. 

R(t,o,o;to,fo,ro)  ♦  A2 (o,to, fQ,ro) 


R(o,n.o;t  ,f  ,rj  J  h*(t;t  ,f  ,r  ) 

°  °  « 

l 

R(o,o,s;to,fo,ro')  *  s2(sin  a;tQ,fo,ro)  . 
sin  a 

where  o,  (,  and  a  are  smear  variables. 

The  extent  of  (nonnegligible)  A2  (the  differential  Doppler  spectrum) 
in  c  represents  the  frequency  smear  imposed  on  a  tone  transmitted  at  f«f0 

by  motions  in  the  propagation  medium. 

Similarly,  the  extent  of  h2  in  C  measures  the  range  of  delays  (time 
smear)  imposed  on  a  signal  due  to  the  different  travel  times  of  the  paths 
connecting  source  and  receiver  (multipath). 

The  extent  of  S2  in  a  is  the  angular  smear  imposed  on  a  transmitted 
plane  wave  by  refractive  and  diffractive  effects  in  the  medium. 

The  functions  A2,  h2,  and  S2  can  be  measured  by  transmitting  signals 
that  are  narrow  compared  with  the  imposed  smears,  corresponding  effectively 
to  impulses  in  time  (a  narrow  pulse),  frequency  (a  tone),  and  angle  (a 
plane  wave) . 

TEMPORAL  ISOLATION  OF  THE  DIRECT  ACOUSTIC  PATH 

.ieasurement  of  the  second  order  statistics  of  the  transfer  function 
R(x,n,s;t  ,f  ,r  )  for  the  direct  acoustic  path  necessitates  isolation  of 
that  path°fr8m  8ther  significant  multipath  arrivals. 

Fortunately,  the  time  occupancy  of  the  direct  path  impulse  response 
was  smal  1  compared  with  its  isolation  in  time  from  nondirect  arrivals. (Thi s 
is  true  for  both  the  18  and  1.5  nmi  paths.)  This  made  it  possible  to 
transmit  pulses  long  enough  to  allow  measurement  of  the  direct 
path  steady  state  Response  but  short  enough  to  prevent  intersymbol 
interference  from  nondirect  paths. 
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Measurement  of  R  also  requires  pulse  repetition  intervals  small 
compared  to  the  direct  path  correlation  time.  It  will  be  shown  below 
that  the  repetition  intervals  used  at  AFAR  (typically  3  and  10s)  yielded 
smooth  estimates  of  the  time  covariance  function.  Although  it  is 
possible  to  contrive  time  varying  transfer  functions  that  will  yield 
smooth,  repeatable  covariance  estimates  even  when  undersampled,  the 
stochastic  nature  of  ocean  volume  propagation  allows  the  conclusion  that 
sample  rates  >0.1  Hi  are  sufficiently  high  at  all  AFAR  acoustic  frequen¬ 
cies. 


EXTRACTION  OF  A  CK  ESTIMATE 

Sequence:*  of  narrowband  Cauaalan  envelope  pulses  were  used  to  sim¬ 
ulate  CW  transmission  over  the  AFAR  direct  acoustic  paths.  The  received 
pulse  Intensities  were  computed  as  time  functions  extending  over  each 
pulse  repetition  interval  and  were  then  averaged  to  determine  the  average 
channel  power  response  over  the  duration  of  each  experiment.  The  trans¬ 
mission  delay  of  the  peak  average  power  response  was  determined  and  the 
Individual  received  waveforms  were  sampled  (tapped)  at  that  fixed  trans¬ 
mission  delay.  The  resulting  sample  sequence  approximates  sampling  a 
CU  transmission  over  an  isolated  direct  path  using  a  uniform  sampling 
Interval  equal  to  the  pulse  repetition  Interval. 

INTENSITY  AND  PHASE  MEASUREMENTS 


The  tapped  sequences  of  quadrature  data  were  used  to  compute  time 
series  of  intensity  and  phase; 


Til 


•  ui  •  >vt 

complex  envelope 

•  u,2  •  V1  2 

intensity 

-  tan*1  Vj/U.  ♦ 

phase 

-  V(2-f0) 

phase  delay. 

ries  of  phase  are  made  continuous  by  finding  such  that 

|#i^l^fNie.l*Ni)  2w  *♦!)  <  e;  N  -0.  (4) 


This  is  Justified  by  assuming  that  phase  variations  are  due  to  smoothly 
varying  pulse  arrival  times  and  (#.  .  -♦,)/2*f  <<  pulse  repetition 

interval.  *  0 
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Figure  3  shows  typical  intensity  and  phase  delay  plots  for  the 
18  nmi  direct  path. 

Empirical  distribution  functions  of  envelope  level  |S^|  were  coa- 
puted.and  typical  results  are  shown  in  figure  4.  At  short  range  and  low 
acoustic  frequency  (weak  scatter  regime), the  statistics  are  log  normal 
as  predicted  by  first  order  perturbation  thoory.11  As  acoustic  frequency 
and  range  are  increased,  the  deterministic  direct  path  begins  to  break  up 
into  a  tube  of  randomly  varying  "micropaths 1  At  sufficiently  high 
frequencies  (and/or  long  ranges),  the  micropaths  add  randomly  at  the 
receiver,  yielding  quadrature  components  with  Gaussian  statistics  (a 
Rayleigh  fading  envelope). 

The  AFAR  1.5  nmi  data  are  apparently  log  normal  and  the  18  nmi 
data  are  Rayleigh*  at  all  center  frequencies. 

A  commonly  used  measure  of  signal  intensity  fluctuation  is  the 
variance  of  log  (I).  First  order  perturbation  theory  predicts  a  log 
variance  that  increases  monotonical ly  with  both  range  and  acoustic 
frequency.1'  This  is  clearly  a  nonphysical  result  and,  in  fact,  the 
variation  in  log  intensity  saturates  at  a  level  corresponding  to  the 
value  predicted  for  Rayleigh  fading.  Figure  5  shows  values  of  log 
variance  computed  using  AFAR  direct  path  data.  The  abscissa  scaling 
is  that  which  will  yield  a  linear  dependence  predicted  by  first  order 
perturbation  theory. Data  In  the  saturated  regime  (Rayleigh  limit) 
were  obtained  on  the  18  nmi  direct  path,  and  data  on  the  Increasing 
portion  of  the  curve  were  obtained  from  the  1.5  nral  path.  The  satura¬ 
tion  curve  is  similar  to  those  observed  In  atmospheric  electronwgnet ic 
propagation  measurements. 


TIME  SMEAR 

The  duration  of  the  time  smear  imposed  on  a  signal  by  a  transmission 
channel  can  be  directly  measured  by  transmitting  pulses  of  short  dura¬ 
tion  compared  to  the  smear.  Figure  6  is  a  set  of  envelopes  of  direct  path 
arrivals  produced  by  1.25  ms  pulses  transmitted  at  a  center  frequency  of  3200  Hr. 
The  multimodal  arrival  structure  is  spread  over  about  a  10  ms  transmis¬ 
sion  delay.  Figure  7  is  an  incoherent  average  of  48000  consecutive 
pulse  arrivals  (a  140-hr  average);  l.e.,  figure  7  Is  an  estimate  of  h^fO. 

A  theoretical  average  power  response  curve  derived  by  the  JASON  group, 
using  measured  oceanographic  parameters.  Is  an  excellent  match.1. 


*  The  JASON  work  predicts  that  the  statistics  are  not  precisely 
Rayleigh  and  investigation  of  higher  order  intensity  moments  has 

verified  this  prediction.11 
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Figure  3a.  Intensity  Tiae  Series 


TIME  IV) 


Figure  3b.  Phase  Delay  Time  Series 

Figure  3.  Intensity  and  Phase  Delay  Tiae  Series,  AFAR  18-nrai  Direct 
Path  (130  as  Gaussian  Envelope  Pulses,  10.158  s  Pulse  Repetition 
Interval),  CW  Approxiaation 
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CUMULATIVE  PROBABILITY 


Figure  4*.  AFAR  Unsaturated  Direct  Path 


CUMULATIVE  PROBABILITY 


Figure  4b.  AFAR  Saturated  Direct  Path 
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Figure  4.  Empirical  Amplitude  Distribution  Functions  for 
AFAR  Unsaturated  and  Saturated  Direct  Paths  Compared  With 
Log  Normal  and  Rayleigh  Distributions  (Plot  Courtesy  H.A.  Freese) 
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Figure  5.  Variance  of  Natural  Logarithm  of  Intensity  Fron 
AFAR  Direct  Path  Data  (Frequency  (f)  in  Hx,  Path  Length  (L)  in  ■) . 
95  Percent  Confidence  Liaits  Are  Shown 
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Figure  6.  Received  Pulse  Envelopes  at  5  »in  Separations 
(1 .25  ms  Gaussian  Envelope  Pulses  at  3200  Hz  Center  Frequency) 
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Figure  7.  h'  (C).  the  Incoherent  Average 

of  140  hr  of  Pulses  Received  Over  the  18-nai  Direct  Path 
(1.2S  bs  Gaussian  Envelope  Pulses,  3200  Hz  Center  Frequency, 
Pulse  Repetition  Interval  of  10.158  s) 
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As  was  pointed  out  earlier,  the  Fourier  transform  of  the  channel 
tine  smear  h*(C;t  ,f  ,r  )  is  the  frequency  covariance  function 
R(o,n,o;t  ,  f  ,r  )?  Pigure  8  shows  the  Fourier  transform  of  lr  (O 

of  figure°7.°  ° 

The  frequency  covariance  function  can  also  be  generated  directly 
by  computing  the  zero  time  lag  crosscovariance  between  samples  of  the 
complex  signal  envelopes,  obtained  at  a  fixed  time  after  transmission, 
of  two  i'W  signals  of  differing  frequency.  The  values  of  covariance  are 
then  plotted  versus  difference  frequency, yielding  the  frequency  covariance 
function.  Figure  9  shows  the  frequency  covariance  functions  computed 
from  several  sets  of  AFAR  direct  path  data  at  various  center  frequencies. 
The  e*1  covariance  widths  are  several  tens  of  hertz  and  compare  favorably 
with  the  independent  measure  of  figure  8.  The  results  show  no  clear 
frequency  dependence, and  the  variation  between  curves  is  probably  due 
to  insufficient  (10  hr)  averaging  time. 


FREQUENCY  (Hf) 


Figure  8.  Frequency  Covariance  Estimate 
(Transform  of  Figure  7,  Time  Spread) 
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Figure  9.  Frequency  Covariance  Estimates  at  Eight  Different  Center 
Frequencies,  AFAR  18-nmi  Direct  Path 


FREQUENCY  SMEAR 

Differential  Doppler  spectra  have  been  computed  for  the  18  nmi 
direct  acoustic  path  bv  transmitting  sequences  of  150  ms  Gaussian 
envelope  pulses  of  fixed  center  frequency.  The  pulse  repetition 
interval  of  10s  is  small  compared  to  the  measured  covariance  times 
at  all  AFAR  experimental  frequencies. 

The  tap  sequences  were  broken  into  50  percent  overlapped,  zero  mean 
segments,  Hanning  windowed,  and  Fourier  transformed.  Magnitude  squared 
spectra  were  averaged  over  all  segments  and  normalized  to  obtain  power 
spectral  density  estimates.  Figure  10  shows  a  typical  result  for  the 
18  nmi  direct  path  obtained  by  averaging  five  40  hr,  50  percent  over¬ 
lapped  segments  (120  data  hr).  Figure  11  contains  smoothed  spectral 
estimates  for  data  taken  simultaneously  at  three  different  center 
frequencies.  The  measured  frequency  smears  are  on  the  order  of  milli- 
hertz  and,  as  expected,  increase  with  center  frequency. 
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The  Fourier  transform  of  channel  frequency  smear  is  the  time  auto¬ 
covariance  function.  Normalized  time  autocovariance  functions  were 
computed  from  the  tap  data  sequences,  using  fast  Fourier  transform 
techniques.  Figure  12  shows  typical  results.  Figure  13  is  a  plot  of 
direct  path  1/e  covariance  time  as  a  function  of  center  frequency 
measured  during  a  number  of  AFAR  experiments.  The  circled  data  arc 
values  obtained  during  the  Joint  Oceanographic  Experiment  of  1975  and 
represent  the  longest  averaging  times  (~100  data  hr).  The  coherence 
time  is  =  1000sat  400  Hz  and  falls  off  approximately  as  f  l. 
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Figure  10.  Differential  Doppler  Estimate,  AFAR  18-nmi 
Direct  Path  (130  ms  Gaussian  Envelope  Pulse,  412  Hi 
Center  Frequency,  10.158  s  Pulse  Repetition  Interval) 
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Figure  11.  Saoothed  Estimates  of  Differential  Doppler  Spread, 
AFAR  18-nmi  Direct  Path  (130  as  Gaussian  Envelope  Pulses, 
10.158  s  Pulse  Repetition  Interval) 


-4000  -2000  0  2000  4000 
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Figure  12.  Autocovariance  of  a  Or  (Approximation) 
Signal  Received  Over  the  AFAR  18-nai  Direct  Path 
(130  ■  Gaussian  Envelope  Pulses,  412  Hx  Center  Frequency, 
110  hr  Average,  10.158  s  Pulse  Repetition  Interval) 
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CENTER  FREQUENCY  |H/) 


Figure  13.  Tiae  Lag  for  Autocovariance  to  Decay  to  1/e  of  Maxima 
as  a  Function  of  Center  Frequencey,  AFAR  18-nai  Direct  Path 
(130  as  Gaussian  Envelope  Pulses,  Pulse  Repetition 
Intervals  <d0.1S8  s).  Solid  Line  is  a  Least  Mean  Square  Fit, 
T1/e  ■  CF~;  C  -  9.2  x  10*.  «  -  -1.12  t  0.0  -  0.06. 


ANGULAR  SMEAR 

Physical  limitations  at  AFAR  Bade  estiaates  of  angular  saear  (or 
its  transfora,  spatial  covariance)  difficult  to  obtain.  Under  certain 
restrictive  assumptions,  it  was  possible  to  derive  estiaates  of  1/e 
covariance  length  for  horizontal  and  vertical  separations  in  a  plane 
perpendicular  to  the  direction  of  propagation.  The  aperture  weighting 
for  each  of  the  antennas  at  the  receiver  end  of  the  18  nai  direct  path 
was  assumed  to  be  Gaussian.  In  addition,  the  spatial  covariance  func¬ 
tion  was  assuaed  to  have  a  Gaussian  shape  in  the  receiving  plane.  The 
assumptions  lead  to  the  following  expression  for  spatial  covariance: 11 
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aB 

C  •  —  oxp 
AB 


C5) 


where 


X  .  Y  are  the  horizontal  and  vertical  separation  of  the  antenna  pair, 

o  o 

a, 8  are  the  i/e  covariance  lengths  in  the  X  and  Y  directions, 
respectively,  and 

“2  -  a2  ♦  p2  r;  ♦  u22K  •  (6) 


'1 


B2  ■  8*  ♦  Rf  ♦  vW2 


(7) 


where  R  is  the  radius  of  the  antenna  (see  figure  1),  and  pR  is  the 
effective  radius  of  the  antenna  aperture  at  which  the  response  has  fallen 
to  i/e  of  the  maximum. 

Since  the  three  antennas  are  siailarly  constructed,  it  was  assumed 
that  p.  ■  p..  Estimates  of  the  horizontal  and  vertical  covariance 
lengths  (a,B)  were  computed  using  the  three  available  antenna  pairs  at 
the  receiver;  the  results  are  plotted  as  a  function  of  center  frequency 
in  figure  14.  The  resulting  covariance  lengths  (»60A  horizontal,  *10A 
vertical)  imply  angular  smears  «1°  horizontal  and  »S°  vertical.  The 
asyaaetry  is  a  manifestation  of  a  corresponding  asymmetry  in  oceanic 
index  of  refraction  variations.*** 
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Figure  14.  Spatial  Separation  for  Covariance  to  Decay 
to  1/e  as  a  Function  of  Pulse  Center  Frequency.  AFAR  18-nni 
Direct  Path  (130  as  Gaussian  Envelope  Pulses,  3  s 
Pulse  Repetition  Interval).  The  Horizontal  Separation 
Is  Shown  With  Dots,  the  Vertical  Separation  With  Triangles 

TIME  LAGGED  SPATIAL  COVARIANCE 

The  covariance  lengths  of  figure  14  were  derived  by  crosscorrelat- 
ing  the  outputs  of  two  antennas  at  zero  tiae  lag.  Since  the  inhoaogenei- 
ties  in  the  ocean  voluae  are  in  aotion,  one  Bight  expect  the  peak  covar¬ 
iance  between  two  antennas  to  occur  at  nonzero  tiae  delay.1 Figure 
IS  shows  tiae  lagged  crosscovariance  functions  coaputed  between  antenna 
P**rs  f°r  the  18-nai  direct  path.  The  (alaost)  horizontally  separated 
antenna  pair  shows  no  evidence  of  a  peak  at  nonzero  tiae  lag,  presumably 
a  reflection  of  the  long  horizontal  covariance  lengths.  In  contrast, 
the  vertically  separated  pair  shows  obvious  nonzero  lag  peaks.  Tie 
dominance  of  vertical  variation  is  consistent  with  an  internal  wave 
viewpoint  as  the  source  of  tiae  variability. **  Figure  16  shows  four 
sets  of  fourteen  6-hr  average  lagged  spatial  covariance  functions  arranged 
in  order  of  increasing  vertical  separation.  Figure  17  shows  the  cor¬ 
responding  result  for  a  single  long  tera  average  (6  X  14  •  84  hr). 
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HORIZONTAL  SEPARATION 
VERTICAL  separation 


HORIZONTAL  SEPARATION  7  6 
VERTICAL  SEPARATION  2  6! 


HORIZONTAL  SEPARATION  7  6 
VERTICAL  SEPARATION  1  15 


Figure  15.  Tine  Lagged  Crosscovariance  Between  Spatially  Separated 
Antenna  (130  as  Gaussian  Envelope  Pulses,  4671  Hi  Center  Frequency, 
10.158  s  Pulse  Repetition  Interval,  6  hr  Average). 

The  Antenna  Pairs  Lie  in  a  Plane  Perpendicular  to  Line  of  Sight. 
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LAGGED  CROSS  SPATIAL  COVARIANCE  FUNCTIONS  8  4671  HZ 


Figure  16.  Tiae  Lagged  Crosscovariance  Between  Vntenna  Pairs 
at  the  Indicated  Spatial  Separations  (130  as  Gauss  it n  Envelope  Pulses 
4671  Hi  Center  Frequency,  10,tS8  s  Pulse  Repetition  Interval). 

AX  Represents  the  Horizontal  Separation,  AY  the  Vertical  Separation. 
Each  Group  Contains  14  Covariance  Estiaates  Froa  Contiguous 
Tiae  Blocks.  Each  Estiaate  is  a  6-hr  Average 
(6  x  14  M  Total  Data  hr). 
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SUM4ARY 

Table  1  suaaarizes  the  aeasureaent  results  presented  for  the  18 
nai  direct  path. 


Table  1.  Stannary  of  Results 


Typical 

Value 

Frequency  Dependence 

Tiae  Covariance  (e'1) 

1000  s  • 

400  Hz 

f-l.l 

Frequency  Covariance  (e "*) 

e 30-50  Hz 

none 

Spatial  Covariance  (e**), 
(Gaussian  shape  assuaptions) 

Horizontal 

>200  a  1 

2  kHz 

f  *; 

«*j  <  «  <  *1 

Vertical 

>4(12 

kHz 

f-H 

In  addition: 

1.  The  envelope  statistics  are  approxiaately  Rayleigh. 

2.  The  variance  of  logf  intensity  is  **/«  (predicted  value 
for  Rayleigh  fading). 

3.  Tiae  lagged  spatial  covariance  functions  reveal  evidence  of 
spatial  transport  due  to  the  notion  of  refractive  index  inhoaogeneities . 

The  results  presented  in  this  report  define  propagation  character¬ 
istics  under  a  very  liaited  set  of  experiaental  conditions.  Their  aain 
values  are  the  stability  of  experiaent  geoaetry,  the  existence  of 
detailed  environaental  data,  and  the  repeatability  of  results.  They 
should  be  useful  in  quantitative  evaluation  of  theoretical  propagation 
aodels.  Clearly,  additional  experiaentation  is  necessary  under  different 
conditions  of  geoaetry  and  environaent.  Finally,  propagation  between 
aoving  platforas  can  be  expected  to  strongly  affect  aeasured  results. 
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EXTRAPOLATION  OF  RESULTS 


This  section  contains  qualitative  extrapolations  of  the  AFAR 
strong  fluctuation  regime  results  to  other  ocean  acoustic  propagation 
conditions.  The  remarks  are  based  on  a  combination  of  measured  ocean¬ 
ographic  parameters  and  physics  reported  by  S.  M.  Platte.13  The  AFAR 
strong  fluctuation  results  presented  above  were  for  a  shallow  angle 
(*  2  deg),  shallow  depth  (3  600  a),  volume  only  path  with  fixed  end¬ 
points.  We  make  the  following  qualitative  general 1  rat  ions  subject  to 
the  caveat,  .  .  acoustic  experiments  Invariably  fall  Into  'anomalous 
regions  .  .  .  . Ml* 


STEEPER  RAYS 

When  the  deterministic  ray  path  departs  significantly  (>  few 
degrees)  from  horizontal,  the  sound  speed  fluctuation  correlation  length 
along/ transverse  to  the  path  dec reases/lncreases  rapidly  due  to  the 
large  ratio  of  horizontal  to  vertical  correlation  lengths  In  the  ocean 
(typically  >  10:1). 

The  contribution  to  acoustic  phase  fluctuation  at  any  point  along 
a  ray  path  decreases  with  the  correlation  length  of  sound  speed  fluctua¬ 
tions  measured  along  the  ray  direction  at  that  point.  Also,  the  number 
of  micropaths  formed  by  scattering  of  a  deterministic  ray  decreases  with 
increasing  correlation  length  transverse  to  the  ray  direction. 

Therefore,  for  a  fixed  path  length,  steeper  rays  have  smaller  rms 
phase  fluctuation  and  form  fewer  micropaths  (l.e.,  longer  paths  are 
required  to  reach  the  region  of  saturated  fluctuations). 


DEEPER  RAYS 

The  variance  of  sound  speed  generally  decreases  with  depth,  result¬ 
ing  In  smaller  acoustic  phase  fluctuation.  Longer  propagation  ranges 
are  required  before  the  phase  differences  In  a  tube  of  rays  along  a 
deep  path  are  large  enough  to  cause  self-interference.  Additionally, 
both  horizontal  and  vertical  sound  speed  fluctuation  correlation  lengths 
Increase  with  depth. 
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GEOGRAPHIC  POSITION 

There  are  certainly  large  variations  in  oceanographic  parameters 
with  geography;  however,  the  limited  experimental  data  available  In  the 
literature  Indicate  some  degree  of  universality  in  spatial  and  tem¬ 
poral  sound  speed  fluctuation  spectra.  Recent  measurements  by  the  AFAR 
group  (unpublished)  reveal  a  similarity  in  the  depth  dependence  of 
sound  speed  variance  for  two  separate  regions  In  the  North  Atlantic 
and  a  third  area  In  the  Culf  of  Mexico. 


PIATFORM  MOTION 

When  either  or  both  endpoints  of  an  acoustic  link  are  In  motion, 
the  correlation  time  of  an  acoustic  path  will  decrease  as  the  translated 
(and  rotated)  path  sweeps  through  independent  regions  of  sound  speed 
fluctuation.  The  large  ratio  of  horizontal  to  vertical  sound  speed  cor¬ 
relation  length  will  result  In  similar  ratio  In  the  relative  effects  of 
horizontal  and  vertical  platform  speeds.  Unpublished  data  of  the  AFAR 
group  show  that  for  platforms  translating  horizontally  and  transverse 
to  an  ocean  volume  acoustic  path,  signal  correlation  time  falls  off 
linearly  with  speed  for  speeds  In  excess  of  ~  5  knots. 


CONCLUSION 

The  information  presented  in  this  document  shows  that  signals 
transmitted  along  acoustic  paths  in  the  ocean  that  contain  no  boundary 
interactions  are  randomly  perturbed  (smeared)  in  time,  frequency,  and 
angle.  Qualitative  extrapolation  to  other  conditions  of  propagation 
is  possible,  and  the  reader  who  is  interested  in  more  physical  insight 
and  numerical  precision  is  urged  to  read  references  11,  12,  15,  and  16. 
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